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Using first-principles supercell calculations, we have investigated energetic, structural and di- 
electric properties of three different A(B'B")03 perovskite alloys: Ba(Zn 1( /3Nb 2 /3)03 (BZN), 
Pb(Zn 1/3 Nb 2/3 )03 (PZN), and Pb(Zr 1/3 Ti 2 /3)03 (PZT). In the homovalent alloy PZT, the en- 
ergetics are found to be mainly driven by atomic relaxations. In the heterovalent alloys BZN and 
PZN, however, electrostatic interactions among B' and B" atoms are found to be very important. 
These electrostatic interactions are responsible for the stabilization of the observed compositional 
long-range order in BZN. On the other hand, cell relaxations and the formation of short Pb-0 
bonds could lead to a destabilization of the same ordered structure in PZN. Finally, comparing the 
dielectric properties of homovalent and heterovalent alloys, the most dramatic difference arises in 
connection with the effective charges of the B' atom. We find that the effective charge of Zr in PZT 
is anomalous, while in BZN and PZN the effective charge of Zn is close to its nominal ionic value. 

PACS numbers: 77.84.Dy, 71.15.Mb, 71.20.Be, 71.23.-k 



I. INTRODUCTION 

Simple perovskite compounds have the chemical for- 
mula ABO3. In the cubic phase, the oxygen atoms form 
a cubic lattice of corner-sharing octahedra with the B 
cations at their centers, while the A cations form a sec- 
ond interpenetrating cubic sublattice located at the 12- 
fold coordinated sites between octahedra. Interestingly, 
most of the perovskite compounds that are of greatest 
technological interests are not simple systems, but rather 
A(B'B")C>3 alloys with two different kinds of B atoms. 
Examples include the Pb(ZrTi)03 alloys that are eptj- 
rently used in piezoelectric transducers and actuatorsou. 
Other examples include Ba(MgNb)C>3 for high-frequency 
applications^, and Pb(MgNb)0 3 and Pb(ZnNb)0 3 which 
exhibit extraordinarily high values af the piezoelectric 
constants when alloyed with PbTiOsH. 

From a chemical point of view, two classes of 
A(B'B")03 compounds can be distinguished: homova- 
lent vs. heterovalent alloys. In homovalent alloys, the 
two B atoms belong to the same column of the periodic 
table. The A(B^.B' 1 '_ :c )03 alloys can thus have a com- 
position x continuously ranging from to 1. A typical 
example of such a homovalent alloy is Pb(Zr x Tii_2;)03. 
On the other hand, there is a unique stoichiometry in 
heterovalent alloys, since the two B atoms do not be- 
long to the same column of the periodic table. Typical 
examples of heterovalent alloys are Pb(Zn!/ 3 Nb 2 /3)03, 
Ba(Zn 1/3 Nb 2/ 3)0 3 , Pb(M gl/3 Nb 2/ 3)03, 
Ba(Mg 1/3 Nb 2/ 3)03, and Pb(Sc 1/2 Ta 1/2 )0 3 . 

While many first-principles calculations have been per- 
formed to understand and predict properties of simple 
ABO3 perovskite systems (see Ref.tTEl and references 
therein) , only a few studies have investigated alloy prop- 
erties via accurate ab-initio techniques. Foe-. example, 
ferroelectric effects in PZT have been studiedc!, and en- 



ergetics of long-range order structures of hamayalent or 
heterovalent alloys have been investigatedlijlij. These 
recent studies provide a useful understanding of the mi- 
croscopic behavior of perovskite alloys, but they do not 
give a simple picture of the consequences of heterova- 
lency on the energetics and properties of the alloys. For 
example, one may wish to compare the driving mecha- 
nisms responsible for the energetics of homovalent and 
heterovalent alloys. It is also unclear how the chemical 
difference between homovalent and heterovalent alloys af- 
fects interatomic distances and dielectric properties. 

Another aspect that is poorly understood is the ef- 
fect of the A-atom identity on properties of heterova- 
lent alloys. For example, many. Ba-compounds in- 
cludin£U3a(Mg 1/3 Nb 2/ 3)03 (BMN)a, Ba(Zn 1/3 Nb 2/ ^CL i 
(BZNP, Ba(Zn 1/3 Tanj3)03 (BZT)BEa, 

Ba(Mg A/3 Ta 2/ 3)03 (BMT)B, Ba(Sr 1/3 Ta 2/ 3)03 
(BST)O and Ba(Ni 1/3 Nb 2/3 )03 (BNN)B, adopt a struc- 
ture with long-range compositional order in which one 
plane of B' atoms alternates with two planes of B" 
atoms along the [111] direction. We shall refer to this 
as the 1:2mh] structure. On the other hand, weak 
X-ray reflections have been detected in various-. Pb- 
compounds, such as Pb(MgnzjNb 2/ 3)03 (PMN)ErEl and 
Pb(Mg 1/3 Ta 2/3 )03 (PMT)E3, that are indicative of a 
rocksalt-like 1:1 ordering of B'-rich and B"-rich sites on 
the B sublattice. We will refer to this as the 1:1m hi struc- 
ture, since it can be regarded as consisting of a simple 
alternation of planes of B'-rich and B"-rich sites along 
[111]. The nature and strength of the atomic ordering 
have been shown to be crucial for the properties of this 
class of alloys. 

Understanding heterovalent and A-atom effects can 
also be very useful for building an accurate effec- 
tive Hamiltonian that can extend the reach of first- 
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principles calculations by peiaining only the physically 
relevant degrees of freedomliil. Effective Hamiltonian ap- 
proaches have been developed in the past for simple per- 
ovskite systems. When constructed from first-principles 
calculations, they have proven to be very successful 
both for reproducing phases transition sequeacesEJcj, 
and for studying ferroelectric domain walla23 as well 
as finite-tesruaerature dielectric and electromechanical 
propertiesclTcJ. An effective Hamiltonian for a per- 
ovskite alloy should include both compositional and 
structural degrees of freedom. A good starting point for 
treating the former may be the simple model of Ref.El 
By focusing on electrostatic effects, this model success- 
fully reproduces the compositional long-range order ob- 
served in a large class of complex heterovalent perovskite 
alloys. However, in order to augment this model with a 
description of structural relaxations, we are motivated to 
understand the role of such relaxations in the energetics 
and the dielectric properties of these materials. 

In the present study, we investigate both heterovalent 
and A-atom effects in an important class of perovskite 
alloys using first-principles techniques applied to ordered 
supercells. To study heterovalent effects, we compare the 
computed properties of Pb(B'B")03 alloys for the cases 
of homovalent vs. heterovalent B-site substitution. With 
regard to A-atom effects, we are mainly interested in the 
differences between Pb compounds and other divalent-A 
compounds. While we find that the energetics of homova- 
lent alloys are mainly driven by atomic relaxations, we 
confirm that the electrostatic interactions are extremely 
important to understand the energetic and other prop- 
erties of heterovalent alloys. For example, electrostatic 
interactions are mainly responsible for the long-range or- 
der seen in heterovalent Ba-compounds. In addition, we 
find that the formation of short Pb-0 bonds in heterova- 
lent Pb(B'B")03 alloys could help in destabilizing the 
structure suggested by electrostatic considerations alone. 
Finally, the present study also emphasizes that the di- 
electric properties of homovalent and heterovalent alloys 
can differ considerably. For example, the effective charge 
of the B' atom is found to be anomalous in homovalent 
alloys, while it is close to its nominal value in heterovalent 
alloys. 

Specifically, we have chosen three compounds for 
study. These are the homovalent Pb(Zr 2; Tii^ :E )03 (PZT) 
system and the heterovalent Pb(Zn!/ 3 Nb2/3)03 (PZN) 
and Ba(Zn 1/3 Nb 2 /3)03 (BZN) systems. To facilitate di- 
rect comparison between the homovalent and heterova- 
lent alloys, we choose the same composition x — 1/3 for 
PZT as for PZN and BZN. Note that in all three alloys 
the larger B atom (Zn or Zr) is the minority species, while 
the smaller one (Nb or Ti) is the majority. 

All the supercells studied in the present article retain 
inversion symmetry, thus preventing the occurrence of 
ferroelectricity and piezoelectricity. While the Pb com- 
pounds PZN and PZT do exhibit ferroelectricity at low 
temperature.— many interesting features (such as short 
Pb-0 bondsHil) appear in the paraelectric phase of these 
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FIG. 1. (a) The l:2 [00 i] A 3 (B'B" 2 )C>9 supercell. The suc- 
cession of planes from bottom to top (i.e., progressing along 
[001]) is (B',0), (Pb,0), (B",0), (Pb,0), (B",0), (Pb,0)... 
(b) B-site atomic ordering in the l:2[m] structure, illustrated 
within a large box aligned with the cartesian axes. The suc- 
cession of planes progressing along the [111] direction is (B'), 
(Pb,0), (B"), (Pb,0), (B"), (Pb,0)...; only B' and B" atoms 
(small open and large filled circles respectively) are shown 
here. 

compounds. Moreover, the atomic ordering on the B sub- 
lattice, which is a focus of the present work, is determined 
at growth temperatures well above the ferroelectric-to- 
paraclcctric transition temperature. For the purpose of 
understanding the intrinsic differences between homova- 
lent and heterovalent alloys, and between heterovalent 
alloys with two different A atoms, we thus decided not 
to relax the inversion-symmetry constraint. 

II. METHODS 

We limit ourselves to the two ordered superlattice 
structures shown in Fig. 1. These 15-atoms supercells 
consist of one B' plane alternating with two B" planes 
along either the [001] direction (Fig. la) or the [111] di- 
rection (Fig. lb). The structure shown in Fig. la will 
be denoted as l:2[ooi]- As previously indicated, l:2[m] 
denotes the structure shown in Fig. lb. The lattice 
vectors of the l:2fooi] structure are arooi] = ao[l,0, 0], 
b[ooi] = Oo[0,l,0], and c [0 oi] = £ a o[0, 0,3], where a Q is 
the lattice parameter and £ is the ordered-related axial 
ratio. The ideal l:2rooi] structure exhibits the tetrago- 
nal 4/mmm point group. On the other hand, the point 
group of the ideal l:2[ m ] structure is 3m, and its lattice 
vectors are ami] = a o[l,0, — 1], b[m] = <Zo[— 1,1,0], and 
c [m] = C a o[l, 1) 1], where ao is the lattice parameter and 
£ is the ordered-related axial ratio. We will retain the 
ideal group symmetry in both structures. The two B" 
atoms are then equivalent by symmetry in both struc- 
tures. Of the three A atoms in the cell, the two A atoms 
located between a B' and a B" plane are equivalent to 
each other and are free to move along the compositional 
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([001] or [111]) direction. These atoms will be denoted 
Ai below. The third A atom, which we denote A2, lies 
on the mirror plane between the two B" planes and is 
forced by symmetry to remain undisplaced. 

We perform local-density approximation (LDA) cal- 
culations on the above supercells p^ing the Vander- 
bilt Jiltrasoft-pseudopotential schemata. As detailed in 
Ref.LJ, a conjugate-gradient technique is used to mini- 
mize the Kohn-Sham energy functional. The ultrasoft- 
pseudopotential approach has two major advantages. 
First, it allows for the generation of exceptionally trans- 
ferable pseudopotentials because of its use of multiple ref- 
erence energies during the construction procedure. Sec- 
ond, a modest plane-wave cutoff can be used, thus also 
allowing us to include the semicore shells in the valence 
for all the metals considered. Technically, we chose the 
plane-wave cutoff to be 25 Ry, which leads to converged 
results of physical properties of interests using roughly 
2,700 plane waves per 15-atom supercell. The Pb 5d, 
Pb 6s, Pb 6p, Zr 4s, Zr ip, Zr id, Zr 5s, Ti 3s, Ti 
3p, Ti 3d, Ti 4s, Zn 3d, Zn 4s, Nb 4s, Nb ip, Nb Ad, 
Nb 5s, O 2s and O 2p electrons are treated as valence 
electrons in the present study. Consequently, our cal- 
culations include 122, 134 and 132 electrons per cell in 
BZN, PZN and PZT, respective!*-. We use the Ceperley- 
Alder exchange and-, correlationt-3 as parameterized by, 
Perdew and Zungerc3. A (6,6,2) Monkhorst-Pack mesht^l 
was used for the l:2[ oi] structure, corresponding to six 
k-points in the irreducible wedge of the Brillouin zone, 
and providing a desirable exact mapping onto the (6J3,6) 
mesh of the primitive 5-atom idxnple perovskite celEj. A 
(3,3,3) Monkhorst-Pack meshc£l was used for the l:2[ ln ] 
structure, corresponding to 10 k-points in the irreducible 
zone of the l:2mi] structure. We checked that this mesh 
leads to converged results by noticing that the total en- 
ergy of PZT l:2[ ln i structure computed with the (4,4,3) 
and the (3,3,3) meshes differed by less than 0.002 eV. 

The lattice parameter ao, the "ordered-related" axial 
ratio £, and those atomic displacements consistent with 
the symmetries of the ideal structures were fully opti- 
mized by minimizing the total energy and the Hellmann- 
Feynman forces, the latter being smaller than 0.04 
eV/A at convergence. 

The dynamical effective charge of an atom was deter- 
mined by calculating the difference in polarization be- 
tween the ground state structure and a structure in which 
the given atom was displaced slightly along the compo- 
sitional direction, the first-order variation of the polar- 
ization with atomic displacement being the dynamical 
effective, charge Z*. We follow the procedure introduced 
in RefEU which consists in directly calculating the spon- 
taneous polarization as a Berry phase of the Bloch states. 
Technically, we use roughly 2,000 Bloch states to assure 
a good convergence of the effective charges. 
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FIG. 2. Illustration of total-energy differences for (a) BZN, 
(b) PZN, and (c) PZT. The zero of total energy is arbitrarily 
chosen to be that of the unrelaxed l:2rooi] structure. 



III. RESULTS: 

A. Energetics 

Figure 2 illustrates the computed energetics for the 
BZN, PZN, and PZT cases, showing the energy differ- 
ences with and without relaxations, and for l:2rooi] ver- 
sus l:2[xn] supercells. The unrelaxed structure has an 
ideal axial ratio £ = 1 and ideal perovskite atomic posi- 
tions, but a lattice constant that minimizes the supercell 
total energy. On the other hand, the lattice constant, 
the axial ratio and the internal coordinates are all op- 
timized in the relaxed calculations. One can clearly see 
that heterovalent and homovalent alloys differ. In both 
of the heterovalent compounds BZN and PZN, the unre- 
laxed l:2mji supercell is lower in energy by more than 
2.5 eV than the unrelaxed 1:2[ 001 ] supercell. In homova- 
lent PZT, on the other hand, neglecting the relaxations 
leads to a nearly identical total energy in the l:2rpoi] and 
l:2r m i supercells. This energetic difference between het- 
erovalent and homovalent alloys can be understood in 
terms of a model recently developed to explain atomic 
ordering in perovskite alloysEl In this model, the total 
energy of a perovskite with B-site compositional freedom 
is approximated by an electrostatic energy given by 



e=— y 



H-l'l 



(1) 



where a is the cubic lattice constant and I runs over the 
ideal positions of the B atoms (e.g., over Zn and Nb atoms 
in BZN and PZN, or Zr and Ti atoms in PZT). e is a 
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electronic dielectric constant providing some screening 
effects and Aqi is the relative charge, with respect to the 
average B charge, of the B atom in cell /. For example, 
since Zr and Ti belong to the same column (IV) of the 
periodic table, Aqi is exactly zero for any cell in PZT. 
Applying Eq. [1] to PZT then leads to an energy that 
is independent of the atomic ordering between Zr and Ti 
atoms. As shown in Fig. 2c, this expectation is nearly ex- 
actly demonstrated in our calculations on the unrelaxed 
l:2[ oi] and l:2[m] structures. On the other hand, in 
both PZN and BZN, Aq is equal to +1 and -2 for Nb 
and Zn atoms respectively. Thus, the electrostatic en- 
ergy E strongly depends on the atomic ordering of the 
B atoms, as demonstrated in Figs. 2a and 2b for the 
unrelaxed hctcrovalcnt supcrcclls. 

One can also evaluate the difference in energy between 
the unrelaxed l:2r 00 i] supercell and the unrelaxed l:2r 11:L ] 
supercell for the heterovalent case (BZN or PZN) using 
the model of Eq. [1]. The result is that the unrelaxed 
l:2[ oi] structure is higher in energy than the unrelaxed 
l:2[m] structure by (5.274 a.u.)/ea. Fitting this for- 
mula to the first-principles energetic difference between 
the l:2[ nl ] and the 1:2[ 001 ] unrelaxed structures leads to 
reasonable values of the dielectric constant e (6.3 in BZN 
and 7.5 in PZN). Thus, analyzing the first-principles re- 
sults in the framework of this electrostatic model sug- 
gests that the screening is larger in Pb-heterovalent than 
in Ba-heterovalent compounds. This explains why the 
energy separation between unrelaxed l:2[m] and l:2[ooi] 
supercells is larger in BZN than in PZN. 

We consider next the relaxation energies (including re- 
laxation of both atomic and lattice degrees of freedom) . 
Is is interesting to note that relaxations provide a com- 
parable decrease in energy (0.6 eV-0.9 eV) in all three 
l:2[ m ] supercells, independently of both the B-atom va- 
lence and the A-atom kind. On the other hand, the re- 
laxation energy in the l:2[ oi] supercell is strongly depen- 
dent on chemical effects. For instance, relaxing atomic 
and lattice degrees of freedom in the Pb-heterovalent al- 
loy PZN leads to a gain in energy as large as 3.2 eV. 
This gain is 8 times larger than the corresponding re- 
laxational decrease in the Pb-homovalent alloy PZT. As 
shown in Fig. 2, the different energetic behavior between 
Pb-homovalent and Pb-heterovalent l:2[ oi] cells has two 
direct consequences: 

(i) The relaxed PZT l:2[ m ] is much more stable (by 
0.4 eV) than the relaxed PZT l:2[ 00 i]- A similar energetic 
decrease of 0.3 eV has also been calculated in Ref.1121 for 
a PZT alloy with 50% of Zr and 50% of Ti, when go- 
ing from l:l[ooi] to l:l[m] B-site cation ordering. (Con- 
sistent with our previous notation, the l:l[ oi] structure 
exhibits a succession of two different B planes along the 
[001] direction.) 

(ii) The relaxation occurring in PZN leads to a destabi- 
lization of the l:2r 11:L i supercell with respect to the l:2[ 00 i] 
one. The gain in energy by atomic and cell optimization 
in the PZN l:2[ oi] structure thus overcomes both the 



electrostatic energy difference of 2.5 eV between unre- 
laxed 1:2[qoi] and 1:2[ U1 ] PZN structures, and the relax- 
ation energy of 0.6 eV in the PZN l:2[ ln ] structureEl 
This may be consistent with the fact that, unlike many 
chemically similar heterovalent alloys, the ground states 
of Pb-heterovalent alloys such as PMN and PMT do not 
experimentally adopt a 1:2 ordering along a [111] direc- 
tion. The destabilization of the l:2[ ni ] supercell has also 
been found in PZNIliJ: relaxed calculations performed on 
five different order structures predicted that the l:2[m] 
supercell was only the third lowest in energy of those 
studied. 

Relaxation in the l:2[ oi] heterovalent BZN supercell 
also lowers the energy by a considerable amount (2.6 eV.) 
However, in contrast with PZN, this large decrease is 
not able to overcome the very large electrostatic differ- 
ence between the unrelaxed l:2[ oi] and unrelaxed l:2[ ln ] 
BZN supercells. This could be explained by our above 
estimates indicating a smaller dielectric constant in BZN 
than in PZN. Our relaxed LDA calculations thus indicate 
that in BZN the l:2[ ln ] supercell is significantly lower in 
energy than the l:2r 001 ], which is consistent with the fact 
that the ground-state structure observed in BZN (and 
in man y similar Ba-heterovalent alloys) is the l:2r 1:L1 i 
onoiffS. 



B. Structural and dielectric effects 

1. A-atom effects in heterovalent alloys: BZN vs PZN 

One may now wonder what are the structural signa- 
tures and the dielectric consequences of the different re- 
laxations that we find in the three materials under study. 
These are summarized in Table I. 

First, Table I shows that the heterovalent alloys ex- 
hibit similar trends in lattice parameters. For example, 
in both BZN and PZN, the axial ratio is larger for the 
1:2[qoi] structure than for the l:2r m ] one. In particu- 
lar, the axial ratio of the 1:2[ 001 ] structure is significantly 
larger than the ideal value of unity (by 5% and 2% in BZN 
and PZN, respectively). This points to the importance 
of lattice relaxation in heterovalent alloys. As usual in 
LDA calculations, the calculated lattice constants agree 
within 1-2% with—the experimental lattice-constants of 
4.096 A for BZNBEl and 4.04 A for PZNEI, the exper- 
imental values having been measured at room tempera- 
ture in the cubic paraelectric phase. 

Furthermore, we find that all the atomic bonds in 
the l:2fixi] structure are very similar in the BZN and 
PZN cases: (i) the A-0 bonds range between 2.80 and 
2.86 A independently of the A atom (i.e., Pb vs. Ba); 

(ii) the Zn-0 bond is lengthened (around 2.05 A); and 

(iii) the Nb-0 bonds are equally divided into a group 
with intermediate lengths (1.91 A) and a group of longer 
bonds with lengths close to that of the long Zn-0 bond. 
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The atomic bonds in the l:2[ oi] structure are also qual- 
itatively similar between PZN and BZN. For example, 
one may notice four different features characterizing re- 
laxed heterovalent l:2[ oi] structures with respect to re- 
laxed heterovalent l:2r m i structures: (i) four short Ai-0 
bonds, ~2.5-2.7 A, that are much-shorter than those ex- 
pected from the sum of ionic radiild (2.70 and 2.82 A for 
Pb-0 and Ba-0 bonds, respectively); (ii) four short Zn- 
O bonds, ~2.0 A; (iii) two very long Zn-0 bonds, ~2.2- 
2.4 A, that are much longer than the 1.95 A expected 
from ionic radii considerations^; and (iv) one short Nb- 
O bond of 1.8 A. 

The differences in atomic bonds between BZN and 
PZN arise mainly in the l:2r oi] structure. For exam- 
ple, some Zn-0 bonds are much longer in BZN 1:2[qoi] 
than in PZN l:2rooi] ; i n particular, the Zn-0 bond in 
BZN 1:2[qoi] is ~0.4 A longer than the value of 1115 
A expected in a purely ionic perovskite compounded. 
Similarly, some Pb-0 bonds in PZN l:2[ oi] are shorter 
by more than 0.1 A than the shortest Ba-0 bonds in 
BZN l:2[ oi]. The formation of very short Pb-0 covalent 
bonds in PZN, and the existence of (unstable) anomalous 
long Zn-0 bonds in BZN, are probably the main quan- 
titative reasons for the larger relaxation energy found in 
PZN l:2[ooi] compared to BZN 1:2 [001] . 

Table I also shows the Born effective charges Z* for 
BZN and PZN supercells, calculated along the compo- 
sitional direction (i.e., along either [001] or [111]). Ctee 
finds similar trends as in simple perovskite compounds^, 
namely values close to +3.0 and +4.0 for Ba and Pb 
atoms, respectively. On the other hand, we find that 
the effective charges of the Nb atoms are much smaller 
(by 1-2 units) in the heterovalent alloy than in the para- 
electric states of simple compounds such as KNDO3 an d 
NaNbOaO. Posternak et al. have shownca that the 
lengths of Nb-0 bonds can have a drastic effect on the 
effective charges: they found that the Z* of Nb decreases 
by more than 2.5 units when going from the cubic para- 
electric phase to the tetragonal ferroelectric phase of 
KNbOs. Thus, the fact that Nb atoms have a much 
smaller effective charge in BZN and PZN than in sim- 
ple paraelectric compounds may be attributable to the 
Nb-0 bond length distribution found in the alloys. Fur- 
thermore, as seen in Table I, the existence of very short 
Pb-0 bonds, which have-been experimentally detected 
in Pb-heterovalent alloysH'Ej, also has a dielectric con- 
sequence: the effective charge of the Pb atom involved 
in those Pb-0 bonds decreases when going from PZN 
l:2[ m ] to PZN 1:2[qoi]- A similar decrease of the Pb 
effective charge has already been seen when going from 
a system with "normal" Pb-0 bonds to a system with 
short Pb-0 bondstl. Similarly, heterovalent l:2[ oi] struc- 
tures exhibit both longer Zn-0 bonds and a decrease of 
the effective charge of the Zn atom. 



2. Heterovalent effects: PZN vs PZT 

As shown in Table I and at variance with heterova- 
lent alloys, the homovalent PZT alloy adopts a near- 
ideal value of the axial ratio in both the l:2[ooi] and 
l:2[m] structures. The relaxation energy occurring in 
PZT l:2fooi] is mainly attributable to the existence of 
longer Zr-0 bonds and shorter Ti-0 bonds, since the 
Pb-0 bonds are very close to their unrelaxed value of 
2.80 A. Unlike the l:2 [0 oi] PZT and the l:2 [m] BZN and 
PZN structures, the l:2[ m j PZT structure exhibits short 
Pb-0 bonds; again, these are presumably responsible for 
the stabilization of the relaxed l:2r m ] PZT supercell over 
its l:2r oi] counterpart. From a dielectric point of view, 
these short Pb-0 bonds lead to an effective charge in 
PZT l:2 [m] that is smaller than in PZT 1:2 [001] . The 
main difference in the effective charges of heterovalent 
vs. homovalent alloys concerns the B' atom (i.e., Zn in 
PZN and BZN vs. Zr in PZT). One can see from Ta- 
ble I that Zr has an anomalously large effective charge 
~6.5) with respect to its purely ionic value of +4. A sim- 
ilar large effective charge of Zr has been found in simple 
compounds and results from the weak hybridization of 
nominally unoccupied Zr 4c? orbitals with the oxygen 2p 
orbitalse3. On the other hand, since the 3d shell in Zn is 
below the Fermi level, hybridization between the O 2p and 
Zn 3d orbitals does not contribute to the effective charge. 
(Incidentally, we find-that this hybridization is strong, as 
reported also in Ref.Ell.) Hybridization between the O 2p 
and Zn Ad states could in principle contribute to Z* , but 
the Zn Ad states are much higher in energy and therefore 
unlikely to contribute strongly. This explains the small 
values of the Zn effective charges displayed in Table I. 

IV. CONCLUSIONS 

We have performed accurate first-principles supercell 
calculations to investigate energetic, structural and di- 
electric effects in heterovalent and homovalent perovskite 
alloys. 

Our main findings are that the energetics of the ho- 
movalent PZT alloy is mainly driven by atomic relax- 
ation and covalent effects. A PZT structure can thus be 
stabilized via the formation of short Pb-0 bonds. The 
formation of these short Pb-0 bonds, which have been 
experimentally seerH, leads to a decrease of the effective 
charge of the Pb atoms. 

On the other hand, we find that the energetics of the 
heterovalent BZN alloys is mainly driven by electrostatic 
interactions among the Zn and Nb atoms. As shown in 
RelEl, these electrostatic interactions lead to the stabi- 
lization of the 1:2 long-range-ordered structure along the 
[111] direction. This 1:2 structure does not exhibit any 
short A-0 bonds, which is consistent with the fact that 
very short Ba-0 bonds have not been detected (to our 
knowledge) in heterovalent perovskite alloys. 
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The situation in PZN is in between those of BZN and 
PZN, in the sense that the optimization of the cell pa- 
rameters and the formation of short Pb-0 bonds in PZN 
could overcome the large electrostatic interaction among 
the Zn and Nb atoms. The destabilization of the 1:2 
structure along the [111] crystallographic direction can 
thus happen in PZN. As in PZT, the formation of these 
short Pb-0 bonds leads to a decrease of the effective 
charge of the Pb atoms in PZN. 

Finally, one drastic difference in dielectric properties 
between heterovalent and homovalcnt alloys concerns the 
effective charge of the larger B-site atom B'. While Zr 
exhibits anomalous large effective charge, the Born effec- 
tive charge of the Zn atom is very close to its nominal 
value. 
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TABLE I. Lattice parameters ao and £ (see text), atomic distances d, and effective charges Z* of Ba(Zn!/ 3 Nb2/3)03 (BZN), 
Pb(Zn!/ 3 Nb2/3)03 (PZN), and Pb(Zr 1; /3Ti2/3)03 (PZT) supercells. Ai is an A atom that relaxes along the compositional 
direction, while A2 is prevented from relaxing from its ideal position by symmetry. B' is the larger B atom (Zr in PZT, Zn in 
BZN and PZN); B" is the smaller one (Ti in PZT, Nb in BZN and PZN.) For each kjeri of bond, the number of them occuring 
per supercell is indicated in brackets. Our calculations agree with the results of Ref.Eil within 0.04 A for all the bond lengths 
of the PZN l:2[m] structure. 



BZN PZN PZT 





L2[ooi] 


L2[iu] 


L2[ooi] 


l:2[in] 


L2[ooi] 


l:2[m] 


a (A) 


4.02 


4.02 


3.97 


4.02 


3.97 


3.97 




1.05 


1.00 


1.02 


0.98 


1.00 


1.00 


d(Ai-O) (A) 


2.66 [4] 


2.85 [3] 


2.54 [4] 


2.80 [3] 


2.74 [4] 


2.60 [3] 




2.91 [4] 


2.85 [3] 


2.87 [4] 


2.82 [3] 


2.81 [4] 


2.81 [6] 




3.28 [4] 


2.86 [6] 


3.28 [4] 


2.85 [6] 


2.91 [4] 


3.07 [3] 


d(A 2 -0) (A) 


2.81 [8] 


2.83 [6] 


2.74 [8] 


2.80 [6] 


2.76 [8] 


2.76 [6] 




2.84 [4] 


2.84 [6] 


2.81 [4] 


2.85 [6] 


2.81 [4] 


2.81 [6] 


d(B'-O) (A) 


2.01 [4] 


2.07 [6] 


1.98 [4] 


2.04 [6] 


1.98 [4] 


2.06 [6] 




2.38 [2] 




2.19 [2] 




2.11 [2] 




d(B"-0) (A) 


1.82 [1] 


1.91 [3] 


1.84 [1] 


1.91 [3] 


1.89 [1] 


1.98 [3] 




2.01 [4] 


2.06 [3] 


1.99 [4] 


2.05 [3] 


1.93 [1] 


2.00 [3] 




2.10 [1] 




2.06 [1] 




1.98 [4] 




Z*(Ai) 


3.10 


3.08 


3.48 


4.44 


4.04 


3.78 


Z*(A 2 ) 


2.49 


2.20 


3.43 


3.34 


3.53 


3.41 


Z*(B') 


2.24 


3.16 


2.59 


2.95 


6.69 


7.76 


Z*(B") 


6.88 


6.54 


7.82 


6.45 


6.65 


5.81 



7 



